Probing many-body states of ultra-cold atoms via noise correlations 
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We propose to utilize density-density correlations in the image of an expanding gas cloud to probe 
complex many body states of trapped ultra-cold atoms. In particular we show how this technique 
can be used to detect superfluidity of fermionic gases and reveal broken spin symmetries in Mott- 
states of atoms in optical lattices. The feasibility of the method is investigated by analysis of the 
relevant signal to noise ratio including experimental imperfections. 



Much of the excitement in the field of Bose-Einstein 
condensation owes to the clear demonstration it provides 
of the wave character of matter. The condensed state of 
bosons involves macroscopic occupation of a delocalized 
single particle state. Consequently, it is characterized by 
sharp density peaks in the freely expanding gas cloud 
after it is released from the trapjlj. Patterns that appear 
when two or more superfluid clouds interfere^, are a 
direct probe of the single particle coherence, amplified 
by macroscopic occupation. 

Recent experiments open intriguing directions for 
studying many body phenomena beyond single particle 
coherence. For example, observation of the superfluid to 
Mott insulator transition as well as experiments in- 
volving ultra-cold fermions near a Feshbach resonance 0, 
address strongly correlated states of matter. The most 
intriguing aspect of such systems is the existence of non- 
trivial correlations and complex order that defy a descrip- 
tion in terms of (single particle) matter waves. Accord- 
ingly, they cannot be characterized simply by the density 
profile of an expanding cloud. For example, the local- 
ized atoms in Mott states of the optical lattice display a 
vanishing interference pattern which can hardly re- 
veal detailed properties of the quantum state. Likewise, 
superfluidity of paired fermions is not evident as a co- 
herence peak in the density profile 5], and detecting the 
order parameter, presents a considerable challenge. Ob- 
servation of some theoretically proposed "exotic" many- 
atom states 0, 0, Q may prove even more elusive. 

In this Letter, we show that the quantum nature of 
strongly correlated states can be revealed by spatial noise 
correlations in the image of the expanding gas. This is 
similar in spirit to measurements of non classical cor- 
relations of light in optical systems 1 1 and temporal 
current noise in mesoscopic conductors |9|- In analogy 
to quantum optics, this technique allows to study mat- 
ter waves that lack single particle coherence. Specifically 
we show that (i) fermionic atoms released from the trap 
would display a clear signature of superfluidty in their 
density correlations. Furthermore, detailed properties of 
the fermionic superfluidity can be determined, such as 
pairing symmetry, and BCS to BEC crossover [ll]: (ii) 
Atoms released from a Mott-insulating state of the op- 
tical lattice display sharp (Bragg) peaks in the density- 
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FIG. 1: Time of flight imaging. The atoms are released from 
the confining potential at time t = and the density of the 
expanding cloud is imaged at a later time. Spatial noise corre- 
lations in this image can be used to probe the quantum state 
at t = 0. 



density correlation function as a consequence of quantum 
statistics; (iii) These peaks can be used to probe the spin 
ordered Mott states proposed for two component bosons 
[3, Il2| . Finally we verify the experimental feasibility of 
the proposed measurements. 

Before proceeding we note that earlier proposals to de- 
tect fermionic superfluidity relied on dynamical response 
of the cloud 01 > inelastic scattering of light to induce and 
measure excitations|l4l llfjj . or to microscopically image 
the pair wave-function in the trap|l5j. In contrast, the 
present technique provides a direct probe of the pair co- 
herence as well as information on the pairing symmetry. 
Earlier work on p robing Mott-insulating phases should 
also be noted |17|. 

We proceed by formulating a detection scheme for 
atoms released from a single macroscopic trap. Suppose, 
for simplicity, that the system is initially in some pure 
state | $ ) . In a typical experimental setup (Fig. , the 
trapping potential is turned off suddenly, and the atoms 
evolve independently under the influence of the free prop- 
agator Uo(t). This is valid provided that the free-atom 
collision cross-section is not too large. Such conditions 
can be achieved by switching the magnetic field to val- 
ues far from the Feshbach resonance when turning off the 
trapQ. 

In such time of flight experiments, the column inte- 
grated density of the expanding cloud is measured by 
light absorption imaging 18] . The images are commonly 
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analyzed by comparing to theoretical predictions for the 
density expectation value: 

(n a (r)) t = ($|f/+(t)^(r)^(r)C/o(t)|$) (1) 

where ip a is the field operator for bosons or fermions and 
a denotes an internal atomic quantum number (spin). 
After a long time of flight the density distribution be- 
comes proportional to the momentum distribution in the 
initial trapped state (n(r)) t (m / ht) (?iQ( r )) . The wave- 
vector Q(r) = mvj (ht) defines a correspondence between 
position in the cloud and momentum in the trap. 

It is important to realize, that in each experimental 
image, a single realization of the density is observed, not 
the expectation value. Eq. is still meaningful, be- 
cause the density is a self averaging quantity. Each bin 
a in the image represents a substantial number of atoms 
N a , while the atomic noise scales as 0(y/N a ). But since 
N a is not macroscopic the density fluctuations are visible. 
They are characterized by the correlation function: 

g a p(v,v') = (h a (r)h (r% - (h a (r)) t (np(r%. (2) 

In analogy with Eq. this can be related to ground 
state momentum correlations: 

G a0 (r,r') OC («Q( r )Q«Q(r')/3> ~ («Q(r)a) (^Q(r')^) ) ( 3 ) 

The proportionality constant is (m/ht) 6 and we shall 
omit it henceforth. We shall be concerned primar- 
ily with pure density-density correlations Q(r,r') = 
E ^^( r ' r ')) which does not require state-selective 
measurement. In practice it may be more convenient 
to consider the quantity An(r,r') = n(r) — n(r') whose 
fluctuations are closely related to Q(r, r'). If (n(r)) t — 
(n(r')) t , then 

(An(r, v'f) t = G(r, r) + G(r', r') - 2<?(r, r'). (4) 

Fermionic superfluids. As a specific example, we now 
consider superfluid states of fermionic atoms. Such su- 
perfluids sustain macroscopic coherence and their trans- 
port properties are similar to their bosonic counterparts. 
However, the first order coherence measured by the den- 
sity profile of the expanding cloud would not reveal their 
superfluid nature. Consider a system at zero temperature 
described by a BCS-like ground state 

\®bcs) =n( Mk + Wka k T a Ui)l°>- ( 5 ) 

k 

Note that this wave function can describe both weak- 
coupling BCS-like pairing, as well as tightly bound 
pairs for which Uk and have a wide momentum 
distribution^]] . The average density profile of the ex- 
panding cloud Q is proportional to the BCS momentum 
distribution function (n(r)) t = 2|«Q( r )| 2 , which is qual- 
itatively indistinguishable from a Fermi distribution at 
T = T C 0. 



The essential difference between these states lies in the 
two particle correlations. For every atom with momen- 
tum k in the BCS state, there is another one at exactly 
— k. This implies, pronounced correlations between den- 
sity fluctuations on diametrically opposite points in the 
expanding cloud. Indeed, a straightforward application 
of (J3J) for the BCS state gives: 

e(r,r') = 2|u Q(r) | 2 | UQ(r , ) | 2 5(r + r'), (6) 

where <5(r + r') is a sharply peaked function of r + r'. This 
sharp peak is a direct analogue of the first order coher- 
ence peak seen in bosonic condensates. It indicates con- 
densation of zero momentum pairs. As for BEC, it has a 
width ~ ht/mL, limited by the finite size L of the initial 
cloud. The pair distribution n s (r) = 2|wQ( r )| 2 |?/Q( r )| 2 
can be read off directly from the weight of the peak. 
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FIG. 2: Fermionic superfluid. (a) Angle integrated weight 
of the correlation peak at diametrically opposite points, (b) 
Weight of the correlation peak for d x n_ y i pairing, (c) and 
(d) (An(r, r') 2 ) at T > (superfluid fraction 0.7). The width 
of the narrow dip ~ M/(mL) is limited by the system's size. 
The dashed line in (d): BEC of tightly bound pairs. 

FfgGI a ) depicts the angular averaged n s (r) (dashed). 
The solid line as well and other plots in Fig [5] depict 
observable column integrated functions. Pairing symme- 
tries other than s-wave, would be detected by higher an- 
gular harmonics in the pair distribution n s (r) and its 
column integral. As an example we plot this function for 
d x 2_ y 2 pairing in Fig. EJb). 

In practice it may be more convenient to measure the 
fluctuations of An(r,r'). Noting that the state (J5J is 
an eigenstate of An(r, — r) with eigenvalue zero, we have 
(An(r, -r) 2 ) = at T = 0. Well away from r + r' = 
(An(r,r') 2 ) - n Q(r) (l-n Q(r) )+nQ (r / ) (l-n Q(r / ) ), which 
is non-zero on the smeared fermi surface of a fermi liquid 
(at T > 0) and for a BCS state. A sharp dip, whose 
width ht/(mL) is limited by the system size, appears 
around r' = r in the superfluid state (Fig[2Ic,d)). Note 
that at finite T the magnitude of the dip is reduced and 
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it disappears at Tc- Heating a tightly bound pair con- 
densate will produce thermal pairs at momenta k > 0, 
producing a wide background dip in analogy to the bi- 
modal distribution in BEC. 

Optical lattice. We now consider atoms initially con- 
fined to an optical lattice. Calculation of (n(r)) and 
(n(r)n(r')) involves operating on the lattice state |<fr) 
with U(tyi/)(r)U(t). In normal ordered expectation val- 
ues this operator can be safely replaced by its projection 
into the lowest Bloch band, A(r, t) = £\ w;(r, t)ai a [l9j . 
Here Wi(r, t) is the free evolution of a Wannier wave func- 
tion, initially centered around the lattice site at Ri and 
a-ia = J drwi(r,0)ip a (r) is a particle annihilation oper- 
ator at this site. At long times Wi(r,t) oc e -lC ^ r '' R< 
times a gaussian envelope, which we shall approximate 
by a square of the same width W = Tvt/(aom). a$ is the 
width of the Wannier state on the lattice. Now (n a (r)) t 
can be related to the off-diagonal correlation function in 
the ground state [lflj: 

(n a (r)) t « (n Q(r) ) = -L £ e^- R ^ {at a a ja }, 

id 

(7) 

The wave- vector Q(r) now defines a correspondence be- 
tween position in the cloud and the quasi-momentum on 
the lattice. Q can also be written in terms of correlations 
in the lattice state: 

g a ,,(r,r') ~ E eiRi< " Q(r)+ ^" Q(r,) < a L^^^^) 

ii'jj' 

+S(r-r')(n(r)) t -(n(r))(n(r')) (8) 

Here S(r— r') is a true (5-function originating from normal 
ordering in the continuum. Note that formulas Q) and 
(JSJ) hold for bosons and fermions. 

In the superfluid state of bosons, where {a\aj) = \^\ 2 , 
n(r,i) exhibits Bragg peaks at Q(r) corresponding to 
reciprocal lattice vectors G. In the Mott state on the 
other hand (ajflj) ~ and there is no interference 
pattern in n(r, t). 

The Mott state, however, displays non-trivial correla- 
tions in the second-order correlation function © associ- 
ated with atom number fluctuations. To illustrate their 
origin, and to make contact with the Hanbury-Brown- 
Twiss effect in quantum optics, we first apply I|7I8|) to 
the "toy" model of two atoms localized in two wells. For 
an initial Fock state | $) = a\a\ | 0) it gives: 

0( r >O = ^cos(-(r-r')) (9) 

where r\ = 1 (—1) for bosons (fermions). The oscilla- 
tions reflect a two body interference effect which amounts 
to bunching in the case of bosons and anti-bunching for 
fermions, as in the textbook 

exampledHIi. 

Addressing the optical lattice we first consider the two 



extremes of a superfluid and a Mott state: 

\SF) =Y[e z <\0), \Mott) =n4=(4) n |0)-(10) 

\/n\ 

According to Eq. ©, Q SF (r,r') = 0, while 

where / is a lattice spacing and ao, the width of a Wan- 
nier function. We conclude that the Mott state sustains 
second order coherence, seen as Bragg peaks in Q(r, r'). 
Note, again, that for a finite lattice of length L, the 
peaks are not true 5-functions. They have a finite hight 
(N/W d f and width £ = 2irht/mL. 

Fig. 01 shows how the two-site bunching signature of 
Eq. 10 develops into coherent Bragg peaks with increas- 
ing lattice size. A fcrmionic insulating state will display 
similar evolution, only the peaks will be replaced by dips. 
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FIG. 3: Normalized density density correlations, g = 
Q{r,r )/{n(r)) t {n(r )) t in the Mott state for a chain of 2,4 
and 40 sites. 

Note that the Bragg peaks do not reflect reduced num- 
ber fluctuations in the Mott state, but rather the quan- 
tum statistics of bosonic or fermionic atoms. In particu- 
lar, we note that a thermal Bose gas on a lattice would 
exhibit similar coherent Bragg peaks, which can be seen 
as a result of the permutation symmetry imposed on the 
wave- function. However in a Mott state with non triv- 
ial internal structure, the noise correlations can reveal 
complex order as we now illustrate. 

As a specific example let us consider spin ordered Mott 
states, that have been predicted in the Mott phase of 
two component bosons on an optical lattice |7J. The 
two components can be represented by a pseudospin-i: 
aj^jj | 0) = 1 1(4) ) • A fairly general mean-field state of 
a spin-ordered insulator is given by the wave-function 

l*q> = II( cos \ a A + ^ sin \ a \l) I 0) • (12) 
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This state has a non zero spin component along the z 
axis and helical component with momentum q on the xy 
plane. Note that the limiting case q = corresponds to 
a xy ferromagnet while q — tt is a Neil state. An ap- 
plication of to the state (|12fl . yields after summation 
over spin indices: 

Q(r, r') oc (cos 4 ~ + sin 4 ~) ]T ~S d (r - r' + -G) 
+ siii 2 ^^i d (r-r'+^H+G)). (13) 

G <r=± ' 

We find that the q ^ spin order produces satellite 
peaks, whose separation from the central peaks at recip- 
rocal lattice vectors is proportional to ±q. The intensity 
ratio of the satellite to the central peaks depends on the z 
polarization of the spins. The satellites reach their max- 
imal relative intensity, half that of the central peaks, for 
spins polarized on the xy plane. Note that in this scheme, 
an x — y ferromagnet cannot be observed directly. How- 
ever a time dependent helical twist can be induced to this 
state by a magnetic field gradient. 

Experimental Issues. To observe the proposed effects, 
two criteria must be met: (i) the atomic noise should 
be observable in an experiment and (ii) correlated fluc- 
tuations should rise above the uncorrelated noise in a 
statistical sense. 

In a typical experiment, © would be averaged over 
narrow cylinders around r and r', whose bases a corre- 
spond to the spatial resolution. The noise of detected 
probe-photons in a bin of area a is a sum of contri- 
butions from atomic and laser light fluctuations. The 
later is fundamentally limited by photon shot noise. The 
atomic noise exceeds the photon shot noise, provided 
the number of photons per bin in the incoming pulse, 
p a > exp(2K)(N (J ) /(r]K 2 ), where k is the optical depth of 
the cloud, which should be chosen close to unity, and r\ 
is the photon detection efficiency The practical limita- 
tion on the number of photons is associated with photon 
recoil which results in image blurring. Under standard 
conditions 10 — 100 photons per atom can be scattered 
in a measurement time of lOus without blurring the im- 
age on the scale of 10/oti indicating that detection 
of atomic noise is possible. A typical cloud of N = 10 7 
e Li atoms j2l| will reach optical depth n ~ 1 after ex- 
pansion to about 1mm. Therefore each 10/im bin will 
contain about 10 3 atoms. The atomic noise contribution 
will cause bin-to-bin variation in the optical absorption 
at the level of few percent, which is detectable with cur- 
rent technology. 

We now verify that in the case of a fcrmionic superfluid 
the correlated fluctuations rise above the background. 
For fixed r the value of An(r,r') 2 will exhibit fluctua- 
tions of order n s (r) as r' varied. To observe the super- 
fluid dip at r + r' = we must average the result over 



a number of independent positions r. Different positions 
yield independent results if they correspond to distinct 
momenta in the trap. Momentum quantization thus lim- 
its the number of independent positions on the 2d image 
of the cloud to (3(A 2 / 3 ), where A is the total number of 
atoms. Consequently the background statistical fluctu- 
ations can be averaged to A -1 / 3 of the dip magnitude 
making it observable in a single image measurement. 

In conclusion, we demonstrated that spatial noise cor- 
relations in the image of an expanding gas cloud released 
from a trap, can reveal key properties of strongly corre- 
lated states of cold atoms. We anticipate that a similar 
technique can be used to observe signatures of more ex- 
otic states such as spin liquids and valence bond solids. 
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